function densidad = armoDensidad(X,Y,Z,sigma)

%Z matriz nroDatos x 2
%X e Y dan la abcisa y la ordenada de los puntos de la grilla con la conv del matu
%sigma varianza de la interp

nroDatos = length(Z(:,1));
densidad = zeros(size(X));

for i = 1:nroDatos
   densidad = densidad + exp(-((Z(i,1) - X).^2 + (Z(i,2) - Y).^2)/(sigma^2));
end

%Normalizacion
pasoX = abs(X(1,1) - X(1,2));
pasoY = abs(Y(1,1) - Y(2,1));

norma = sum(sum(densidad))*pasoX*pasoY;
densidad = densidad/norma;
